STABILITY ESTIMATES FOR THE RADON TRANSFORM 
WITH RESTRICTED DATA AND APPLICATIONS 



PEDRO CARO, DAVID DOS SANTOS FERREIRA, AND ALBERTO RUIZ 

Abstract. In this article, we prove a stability estimate going from the 
Radon transform of a function with limited angle- distance data to the 
L p norm of the function itself, under some conditions on the support of 
the function. We apply this theorem to obtain stability estimates for an 
inverse boundary value problem with partial data. 
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1. Introduction 

In this work we prove a stability estimate from the Radon transform 
with limited angle-distance data to a local L p -norm of the function. Our 
original motivation to study this problem was to obtain stability estimates 
for the inverse problem in electric impedance tomography (E.I.T.) proposed 
by Calderon. Nevertheless, we think that the results obtained on the Radon 
transform restricted to some partial data sets are interesting by themselves 
and are the main contribution of this work. 

Calderon's inverse problem deals with the recovery of a conductivity 7 in 
the interior of a smooth domain $7 from boundary measurements realized 
by the Dirichlet-to-Neumann map. Let u be the solution of the Dirichlet 
boundary value problem 

Jdiv^Vtt) = in n 

[u\aa = feH^(dn) 

where 7 is a positive function of class C 2 on Q. The Dirichlet-to-Neumann 
map assigns to a function / £ H 2 (<9f2) on the boundary the corresponding 
Neumann data of (jl.ip 

A 7 / = jd u u\an 



(1.1) 
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where d u denotes the exterior normal derivative of u. This is a bounded 
operator A 7 : Hz(d£l) — > H~?(d£l) — in fact a pseudodifferential operator 
of order 1 when 7 is smooth. The inverse problem formulated by Calderon 
[12] is whether it is possible to determine 7 from A~. In fact in its initial 
formulation, the problem concerns only positive measurable conductivities 
bounded from above, and it was solved in dimension 2 in this degree of 
generality by Astala and Paivarinta [I] and remains so far open in higher 
dimensions. 

This question is related to the inverse problem of determining a bounded 
potential q £ L°°(0,) in the Schrodinger equation 



from boundary measurements. This reduction was exploited by Sylvester 
and Uhlmann in [43] and in combination with the boundary determination 
results on the conductivity obtained by Kohn and Vogelius [M] allowed 
them to solve the Calderon problem for smooth conductivities in dimension 
n > 3. When is not a Dirichlet eigenvalue of the Schrodinger operator 
— A + q, the measurements are implemented by the Dirichlet-to-Neumann 
map, which can be similarly defined as for the conductivity equation by 



With a slight abuse of notations, we use the convention that whenever 
the subscript contains the letter q, the notation refers to the Dirichlet-to- 
Neumann map related to the Schrodinger equation (jl.2p , while if it contains 
the letter 7 it refers to the map related to the conductivity equation (jl.ip . 

In the inverse problem with partial data one wonders whether one or the 
other of the Dirichlet-to-Neumann maps A 7 ,A g measured only on a subset 
of the boundary, determines the conductivity 7 or the electric potential q 
inside ft. In dimension two, this problem is settled by the articles \29\ [30] 
of Imanuvilov, Uhlmann and Yamamoto using ideas from Bukhgeim |10] 
who dealt with the inverse problem for the Schrodinger equation with full 
data. See [23] for this problem on Riemann surfaces. In dimension higher 
than three, the first results were obtained by Bukhgeim and Uhlmann [11] 
but required measurements on roughly half of the boundary. The results 
obtained by Kenig, Sjostrand and Uhlmann [33] are the most precise so far 
in dimension n > 3 since they require measurements on small subsets of the 
boundary for, say, strictly convex domains 0. This result has been extended 
to the Dirac system by Salo and Tzou in [39J. We should also mention the 
local inverse problem, in which all the measurement are restricted to input 
Dirichlet data supported on the same (the accesible boundary) subset as the 
output measurements. This problem was settled by Imanuvilov, Uhlmann 
and Yamamoto |29[ [30] in dimension n = 2 and only for very special cases 
(the complement of the accesible boundary being a piece of a plane or a 
sphere) in dimension n > 3 by Isakov [31] and extended to Maxwell equation 



(1.2) 




A g / = d u u\ 9 a- 
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in |13] and |15j . The linearized inverse Calderon problem with partial data 
was studied in [20] . 

Let us describe Bukhgeim and Uhlmann result in more details. For this 
purpose, given a direction £ G S n_1 , we consider the ^-illuminated face of 

an 

= {x e 90 : ($,i/(a:)) < 0} 

and the ^-shadowed face 

dn + (£) = {xedfl: (£,v(x)) > 0}, 
where v(x) is the exterior normal vector at x. 

Theorem 1.1 (Bukgheim and Uhlmann Let £1 be a bounded open 

set in R n ,n > 3 with smooth boundary and let us consider F C dCl an 
open neighborhood of the face dCl-(£). Let qi,q% be two bounded potentials 
on Cl, suppose that is neither a Dirichlet eigenvalue of the Schrddinger 
operator —A + q\ nor of —A + q<i, and that for all f G H l / 2 (dCi) the two 
Dirichlet-to-Neumann maps coincide on F 

= A 92 /|ir, 

then the two potentials agree q\ = q%- 

To describe the uniqueness result of Kenig, Sjostrand and Uhlmann we 
need to introduce the appropriate parts of dfl. Assume yo is n °t in the 
convex hull of Cl, we define the yo-ill umina ted face as 

dfl^(y ) ={i£ dCl : (x — y , v{x)) < 0} 

and the yo-shadowed face as 

dn + (y Q ) = {x G dCl : (x — y , v{x)) > 0}. 

Note the abuse of notation when writing dd±(£) and dCl±(yo), since the for- 
mer one denotes the ^-illuminated and ^-shadowed faces from the direction 
£ while the latter one denotes ^-illuminated and ^-shadowed faces from the 
point yo- Then 

Theorem 1.2 (Kenig, Sjostrand and Uhlmann [33]). Let CI be a bounded 
open set in R n ,n > 3 with smooth boundary and let us consider F,B C dCl 
two open neighborhoods respectively of the faces dd-(yo) anddCl + (yo). Let 
q\,q<i be two bounded potentials on Q, suppose that is neither a Dirichlet 
eigenvalue of the Schrddinger operator — A + q\ nor of — A + qi, and that 
for all f G H l l 2 (dQ) supported in B the two Dirichlet-to-Neumann maps 
coincide on F 

A 9i /|f = Ag 2 /|jr, 

then q\ = q2- 

The main goal of this article is to derive an estimate for the Radon trans- 
form which yields the corresponding stability estimates for the above unique- 
ness results (actually only a local estimate in the case of Theorem 1.2). This 



4 



estimate, which we call a quantitative version of Helgason-Holmgren theo- 
rem, will be obtaned in section 2 (see Theorem 12.51) . Stability estimates for 
the conductivity inverse problem in dimension higher than three go back to 
Alessandrini's article [TJ. This was followed by results in two dimensions by 
Liu [35], Barcelo, Barcelo and Ruiz [5], Barcelo, Faraco and Ruiz [6] and 
finally by Clop, Faraco and Ruiz [18] for discontinuous conductivities cor- 
responding to the uniqueness results of Astala and Paivarinta [4J (see also 
[21j). Other stability results for the Calderon problem in dimension greater 
than two are [25] and |16j . In the case of Maxwell equations the stability 
was obtained in [14j . 

Concerning the inverse problem with partial data, stability estimates cor- 
responding to the results of Bukhgeim and Uhlmann were derived by Heck 
and Wang [23], and in the presence of a magnetic field by Tzou [33] . We 
mention also the uniqueness results obtained by Ammari and Uhlmann [3] 
in the case where the potential is known close to the boundary, and the 
corresponding stability estimates obtained by Fathallah [22] and Ben Joud 
[7J . One single log stability estimate was obtained by Alessandrini and Kim 
[2J in the case of the conductivity equation when the conductivities coincide 
on a neighborhood of the boundary with a known one. The stability of 
the local problem under similar condition as in [31] was proved by Caro for 
Maxwell equations in |15j . 

Let F and B be boundary neighborhoods of the illuminated and shadowed 
faces respectively. The natural norm to consider on the partial Dirichlet-to- 
Neumann map is 

\\K q \\ B ^F = sup I (k q (<p) V>) : IM|jrV2(en) = = 

suppy? C B,suppifj C Fj, 

where (-|-) denotes the duality between H 1 / 2 ^^}) and H~ 1 / 2 {d^L). We'll also 
have to consider a larger norm related to solutions of Schrodinger equation 
belonging to the space H{$1, A) (see section 3.1). This norm was considered 
by Nachman and Street in [37] , where they prove the reconstruction of some 
2-plane integrals of the potential from partial data. We will denote this 
norm as 

The class of allowable potentials under consideration will be in Besov spaces 

K(M,\,p) = {q£ L°°(n) and ql n G W X > p (R n ) : + \\q\\ w x, P < M}, 

where A > 0. This class of potential has the advantage of allowing very 
rough functions if A is sufficiently small. Our stability results are as follows. 

Theorem 1.3. Let £1 be a bounded open set in R n , n > 3 with smooth bound- 
ary. Given an open set N in S™ -1 consider F,Bc dQ, two open subsets of 
the boundary which are respective neighbourhoods of the faces <9f2_(£) and 
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<9f2_l_(£) for all directions £ G iV. Given M > i/tere exists a constant C > 
swc/i t/iat t/ie following estimate holds true 

(1.3) lift - g 2 || LP < C (log | log ||A 5l - A q2 \\* B _ F \y X/2 

for all allowable potentials gi,(/2 G JC(M, X,p) on Q, with 1 < p < oo and 
< A < for which is neither a Dirichlet eigenvalue of the Schrddinger 
operator —A + q± nor of —A + q^. 

Next we consider the case of illumination from a point. Let N be an 
open set which does not cut the convex hull of f2. We will define P, the 
convex penumbra boundary from N, as the set of points x G dQ such that 
there exist a y G N with {x — y,v(x)) = and the hyperplane through x 
normal to v(x) being a supporting hyperplane of ft. In order to keep the 
exposition simple, and relate to the Radon transform (rather than the two- 
plane transform in high dimensions) we will restrict ourselves to the three 
dimensional case n = 3. 

Theorem 1.4. Let Q be a bounded open set in R 3 with smooth boundary. 
Given an open set N in R 3 which does not cut the convex hull ofQ, consider 
two open subsets F, B of the the boundary which are respective neighbour- 
hoods of the faces dfl-(y) and dQ + (y) for all y G N. Given M > 0, there 
exist an open neighborhood G C R 3 of the convex penumbra P and a constant 
C > such that the following estimate holds true 

(1.4) ||gi - q 2 \\LP(G) < C (log | log II V - A g2 || B ^ F |) 

for all allowable potentials q\,q2 G /C(M, X,p) on £1, with 1 < p < oo and 
< A < 1/p, for which is neither a Dirichlet eigenvalue of the Schrddinger 
operator —A + q\ nor of —A + qi- 

The proofs of these theorems will be carried out by using the approach 
of [19], which uses the Radon transform. One can see that the result on the 
Radon transform is general enough to be applied to get stability for partial 
data in the context of [33] in dimension three (the Dirichlet-to-Neumann 
map in this case controls the 2-plane transform, which is indeed the Radon 
transform in three dimensions). This can be achieved with the natural 
norm H l l 2 {dVt) H~ 1 / 2 (dQ,) of the Dirichlet-to-Neumann map, by using 
the solutions constructed by Chung in [T7] . 

Theorem 11.31 was proved in [23] by Heck and Wang, without the condition 
of the Dirichlet data being supported on B and the norm in the partial 
Dirichlet-to-Neumann map considered from H 3j/2 (d£l) to H 1 / 2 ^^!) instead 
of the norm || • 11^^.^. They use the Fourier transform. The change to the 
Radon transform illustrates the use of Theorem 12.51 

The structure of this paper is as follows. In section [21 we will state 
and prove the theorem for the Radon transform which is the main result 
in this work. Section [3] and Section [4] are devoted to prove Theorem 11.31 
and Theorem 11.41 applying the stability estimates for the Radon transform 
proven in Section [2j 
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2. Stability for the local Radon transform 

We recall that the Radon transform of a continuous compactly supported 
function / is given by 

Kf(s,uj) = J S((x,cj) - s)f(x)dx, s G R, u G S n ~ l . 

It is always possible to define, by duality, the Radon transform on compactly 
supported distributions. For the time being, we content ourselves with con- 
tinuous functions with some kind of decay, but later on we will extend its 
definition to a wider class of functions. This transform is even. It is some- 
times convenient to think of the Radon transform as a function on II n_1 , 
the Grassmannian set of hyperplanes in R n . Let H belong to II n_1 , then 

nf(H)= [ f(x)dfi H (x) 

where dfin is the Lebesgue measure on H. To relate both notations in a 
coherent way, we set 

H (s,cj) = {x G R n : (x,lj) = s}. 

Note that w is a unit normal to the hyperplane Hq(s,co) and s a signed 
distance to the origin. For later convenience we might change the origin of 
the affine reference for the Radon transform to the point yo £ R n - If ° ne 
describes H as 

(2.1) H = H yo {s,co) = {i£ R n : (x-y ,oj) = s} 

for some ui G S n_1 and s G R. Relating u, s to H as above, one can define 

(2.2) K yo f(s,u)= [ fdfi H . 

We will also make use of the following notation 

H±(s,u) = {x G R n : ±((x - y ,oj) - s) < 0} 
to denote the half-spaces delimited by H yo (s,uj). 
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We refer to Helgason's book [26] for a general study of the Radon trans- 

form. Of particular importance is the issue of local inversion of the Radon 
transform: given a function / with some a priori regularity and some decay 
at infinity, such that TZf(H) = for every H S S C II n_1 , does / vanish 
on E = Uh£sH? For instance, the celebrated Helgason's support theorem 
reads as follows. 

Helgason's support theorem. Let f be a rapidly decreasing continuous 
function such that its Radon transform vanishes on all hyperplanes disjoint 
from a compact convex set K 



then the support of f is contained in K. 

We are interested in the microlocal approach (which differs from Hel- 
gason's original proof) to prove Helgason's support theorem presented in 
[U [U [28] . This approach is somewhat flexible since it does not require the 
full family of hyperplanes used in Helgason's theorem but can be adapted 
to provide weaker support results when the Radon transform only vanishes 
in a neighbourhood of a fixed hyperplane. A result that follows from this 
approach is : 

Microlocal Helgason-Holmgren Theorem. Let f be a compactly sup- 
ported continuous function such that its Radon transform vanishes in a 
neighborhood of H ((x , £ ), Co)- If supp/ C flo~((xo,£o),fo) then x £ 
supp/. 

From the inverse problems point of view, the above result was used in [19] 
to prove the unique determination of the electric potential and the magnetic 
field in a magnetic Schrodinger equation from partial data. It served as a 
substitute to the original but somewhat more involved argument of Kenig, 
Sjostrand and Uhlmann in [33] also based on analytic microlocal theory. 
Similar ideas were used in [20] to investigate a linearization of the Calderon 
problem with partial data. 

The main result in this section, Theorem 12.51 is a quantitative version of 
the microlocal Helgason-Holmgren theorem. We want to relax the compact 
support and continuity assumptions on /, in order to apply the correspond- 
ing results to the study of the stability of Calderon's inverse problem. In 
the next paragraphs we review some concepts and results that will be basic 
in the microlocal approach and that will clarify the proof of Theorem 12.51 

2.1. Microlocal Helgason's support and Kashiwara's Watermelon 
theorems. We will use the classical notation w 2 = w\ + • • • + to denote 
the holomorphic continuation of the Euclidean scalar product — particularly 
to avoid confusion with the norm \w\ 2 = |^i| 2 +- • -+|7/; n | 2 of complex vectors. 
The Segal-Bargmann transform of an L°° function is given by 



nf(H) = o, hhk = 




c 
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Note that it has the following exponential growth 

(2.3) \%f(z)\ <(27r/ l )fe^( Imz ) 2 ||/|| L oc. 

By duality, it is easy to extend this transform to tempered distributions. 
This transform has a wide range of applications in Analysis; amongst others, 
it provides a way of describing analytic singularities of a distribution^ on an 
open set Q C R n . 

Definition 2.1. A distribution f G T>'{0) is said to be microlocally expo- 
nentially small at (xo,£o) & T*Q if there exist a cutoff function \ & Co°(^) 
such that x( x o) 0; ^ wo constants c, C > and a neighbourhood V ZQ of 
zq = xq — i£o * n C n such that the following improved bound holds on the 
Segal-Bargmann transform 

(2.4) \T h ( X f)(z)\<Ce-f + ^ lm ^ 2 . 

for all z £ V ZQ and all h £ (0, 1]. 

The analytic microsupport of a distribution f - - which we denote by 
M- su PPa / — *s ^ e complement of the set of covectors (xo,£o) £ T*£l at 
which f is microlocally exponentially small. The analytic wave front set 
WFa / of f is the complement in T*Q \ of the set of covectors at which f 
is microlocally exponentially small. 

The analytic microsupport is a closed conic0 set of the cotangent bundle 
and consists in two parts 

^supp A / = supp / x {0} U WF A /. 

The projection with respect to the space variable of the analytic wave front 
set is the analytic singular support 

vr(WF A /) = suppa /, vr : T*Q -> fi 

i.e. the set of points xq £ R n which have no neighbourhood on which / 
is real analytic. A microlocal form of Helgason's support theorem reads as 
follows. 

Microlocal Helgason's theorem. If the Radon transform lZf(s,oj) of 
f € C°(R n ) vanishes in the neighbourhood of (sq,loq) G R x S™ -1 then 
(x ,oj ) £ WF A / where x G H (s ,ujo). 

A more invariant formulation would be that the conormal N* (Hq(so, ujq)) 
to the hyperplane is contained in the complement of the analytic wave front 
set. As we will see it implies a local weaker (but quite flexible) form of 
Helgason's support theorem. 

In the situation where a distribution is supported on one side of a hyper- 
plane, Kashiwara's Watermelon theorem describes some of the covectors of 

^In fact, this analysis can be extended to hyperfunctions. 

2 That is, in the frequency variable: (x,£) € u.supp A / =>• (x, A£) £ (i.supp A / for posi- 
tive A. 
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the analytic micro-support ([32], [JTJ Theorem 8.3.3], [271 Theorem 9.6.6], 

EH). 



Kashiwara's Watermelon theorem. Let f G D'(R n ) be a distribution 
supported on one side H + of a hyperplane H . Let uq denote a unit normal 
to H . If (xq, £o) £ M-suppa / then so does (xq, £o + tvo) for all f 6 R. 

Kashiwara's Watermelon theorem is generally stated in terms of the an- 
alytic wave front set, we chose to use a formulation involving the analytic 
microsupport (for similar formulations, see also [5B]) since it encompasses 
information about the support of the function. In particular, it immediatly 
implies the following unique continuation property 

supp/ C {{x -xq,vq) < 0} and (x ,u ) g WF A (/) =4> x £ supp/. 

This is sometimes known as Holmgren's microlocal uniqueness theorem^ (or 
the co-Holmgren theorem). 

2.2. Relating the Radon and the Segal-Bargman transforms. In this 
paragraph, we want to connect the Radon and the Segal-Bargman trans- 
forms. We start from the identity 

/(era;) = Kf(a,uj) 

where lZf(a, co) denotes the (one-dimensional) Fourier transform of 1Zf(s, uj) 
with respect to s 

/oo 
e- isa Kf(s,oj)ds 
-oo 

and use PlancherePs identity to compute the scalar product 

/ fgdx = (2vr)" n / / TZf-(a,uj)1Zg(a,uj) a 11 ^ 1 da du. 
J Jo Js™- 1 

Using the fact that the Radon transform is even, and once again Plancherel's 

identity, we get 

(2.5) [fgdx = \{2*r n+1 f°° [ Kf(s,co)\D\^Kg(;oj)(s) dsduj. 

We choose g to be the conjugate of the Gaussian kernel of the Segal-Bargman 
transform: we begin by computing its Radon transform 

TZ(e-^- x)2 )(s,oj) = f S((x,u) -s)e-^ 2 - x ^dx 



n-l 



(2irh) 2 e 2h 



£(*-<".*» 2 



In the usual formulation using the analytic wave front set, this is obtained by a limiting 
argument from the Watermelon theorem. If xq G supp / then / cannot be analytic at xo 
hence there exists (a;o,£o) G WFa / therefore (xo,£o + tvo) G WFa /■ By the conicity 
of the wave front set, we have (xo,£o/t + uo) G WFa / and by the closedness, we have 
(x ,vo) G WF A /. 
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and plug g = e~2E ( - z ~ x ' in the identity (|2.5p to compute the Segal-Bargman 
transform of a function in terms of the Radon transform 

1 -i-i f°° f 

(2.6) T h f(z) = -(27r)-^-h^- / G n (8,(u,z))Hf(8,u)d8du 

where the kernel G n is given by 

G n (s,w) = \D\ n - 1 (e-^- w ^)(s), sGR, w e C. 
We will use the following estimates of the kernel: 
Lemma 2.2. The kernel G n satisfies the following bound 

(2.7) |G n (s,u>)| < 3^-^(1 + h-^\s-w\) n (l + e^ Imw ^^ s - Rew ^). 

Proof. We need to distinguish two cases according to the parity of the di- 
mension n. 

Let us start with the case n odd which is simpler. The kernel G n can 
explicitly be computed 

G n (s,w) = D^-^e-^ 3 -^ 2 ) = bT^e-^'-^Q. 

where 

Q n (w) = e^D^e- 3 ^ 
is a Hermite polynomial of degree n — 1, hence satisfies the bound 

\Qn(w)\ < A n (l + M) n-1 . 

The former estimate together with 

\ e -jK(s-™) 2 \ = e -27r( s - Rc -) 2 +27r( Im -) 2 
imply the following bound on G n 

(2.8) \G n (s, w)\ < A n h-^ (1 + br*\a - w \) n - x 

when n is odd. 

Notice that identity (|2.6|) reads in odd dimensions 

T h f(z) = 

i(2vr)"^ l^l^e-^^^Q^^^^/^^dc.d,. 

The even dimensional case is a bit more involved: the kernel G n satisfies 
the following relations 

(2.9) G n (s,w) = G n (s — Rew,ilmw), 
G n (s,w) = G n (s,w) = G n (-s, -w) 
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and has the following expression 

G n (s,w) = \D\ n - 1 (e-^- w ^)(s) 

^e"^ [°° _! _ (^/^)\ i s 
= h 2 / \a\ e 2 V s dcr. 

V27T J_oo 

By the relations (|2.9p . it suffices to prove the estimate when w = —it is 
imaginary and s is non-negative, and by scaling, we might as well assume 
that h = 1. For r G R, s > and n even, we may decompose the kernel as 

r 2 ■ 

e ' ' / I , n _l _£__l_i scr 



G w (a,-ir)= ( / -/ |cT-r| n - i e-- +4S,T d ( T 



= /„(s,r) + I n (s, -t) 

where the integral 



r 2 

1ST PT 2 

I n (s, r) = / (r - ^""ie-V-H^ da 

V27T J_oo 

can be computed by an integration on the contour (— oo, r] U [r, r + is] U [r ■ 

is, — oo + is) 

e 2 / f T , (a + is) 2 , ■ , | - \ 

I n (a,r) = -^=( / ( T - a - zs) n - 1 e-^-^ +ls ( CT+IS ) da 



+ (-z) n jT fJ ™"i e -5( T + i<T ) 2 e JS ( T + i<T ) da^j . 

The first term is bounded by a constant times 

e 4-4(l + \ s \ + \T\) n ~ 1 
while the second term3 is bounded by a constant times 

(i + N) n 

This completes the proof of the lemma. □ 
We restate (|2.6[) in term of the Radon transform centered at yo, see (|2.2p . 



as 



(2.10) 75i/(C) 
l 

2 



-(2tt) 'V/t'V / / G„(s, (w,C-yo))^ /( s > w ) dsdw - 



oo JS 



4 In the odd dimensional case, this term disappears when one computes J„(s,r) + 
I n (s, -t). 
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Given ojq E S n and (3 G (0, 1] we consider the following cap centered 
around ujq on the hypersphere S n_1 



(2.11, 

' - on— 1 



r = {a;GS rt - 1 :(w,w ) 2 >l-/3 2 } 
= {u> E S n_1 : d S n~i(ujQ,uj) < arcsin/3} 



d S n-i being the geodesic distance on S n_1 . 

Before proceeding to further computations, we also note that d[j,H (s,Lj) A 
ds = dx and therefore 



/oo 
\Kf(s,u)\ds < 
-oo 



L 1 



L 1 • 



which leads to 

/oo /• 
/ dwds < IS™" 1 ] > 

-oo Js™- 1 

We introduce the following set of functions: u £ X by definition if and 
only if u € L 1 (R n ) and 

(2.13) \Hx= [ (l + |s|) n ||^ n(s,-)|| L i { s-i) ds<oo. 

J R, 

Let us remark, see (|2.12p . that for functions in L 1 (R n ), the Radon trans- 
form is defined a.e. as a function in L X (R x S n_1 ). Our space X is more 
restrictive, a sufficient condition for a function to be in X, is given by the 
estimate 

(2.14) \\u\\ x < IS™" 1 ] / (l + \x\) n \u(x)\dx. 

Proposition 2.3 (Quantitative Microlocal Helgason's theorem). Let f be- 
long to X . There exists a positive constant C, only depending on n, such 
that 

(2.15) e-^ l ^ 2 \T h f(C)\ < £(1 + |C| + \yo\T 

h'2 

[ (1 + \s\) n \\K yo f(s, -)|| Ll(r) ds + \\f\\ x fe-i + e"^) Y 

J \s\<Ot / 

for all h E (0, 1], 7 > and ( G C n such that \Re( - y \ < a/2, |Im(| > 7 
and ( W0 ,^) 2 >l-/3 2 /4. 

Remark 2.4. Note that when TZf vanishes on a neighborhood of (sq,ujq) 
the above estimate implies that (yo,±uJo) £ WFa(/) when (y ,wo) = sq. 
This is the microlocal version of Helgason's support theorem as stated in 
the introduction of this section. Proposition 12.31 is therefore a quantitative 
version of this microlocal result. 
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Proof. It follows from Lemma 12.21 and (|2.10p that 
\T h f(()\<C [°° [ (l + ^l| s -( w ,C-yo)|) n 

J-oo JS™- 1 

x (i + e 27r((^im0 2 -(^-KRcC-s/o)) 2 )) \n y J(s,uj)\duds. 
Let us split the integral into 



+ / + / 

'RxS"- 1 7R\(-d,«)xS"- 1 J(-a,a)x(S™- 1 \r) J(-a,a)xT 

To estimate Ii, notice that if s G R \ (—a, a) and ui G S n ~ , then 

e _^ Ka , iReC _ yo) _ s |2 < e -^j(|.|-|( W ,ReC-tt)>|) a < e -27r4, 

and that 

l< e -^ 2 e2kl Im ^l 2 , 

for all C £ C n such that |Re£ — yo| < a/2 and |Im£| > 7. Then it follows 
easily that 

h < ^e^l Im <l 2 (l + |C|) n ll/IU + e"rf). 

n 2 

The integral I2 can be bounded in a similar way. Notice that, for s G (—a, a) 
and w G S™^ 1 \ T, it holds with 6 = (ImCr^mC 

e5 L(^ImC> 2 < e ^|ImC| 2 e -^(|ImC| 2 -(^ImC) 2 ) < e ^|ImC| 2 e -^7 2 (l-<^} 2 ) ; 

and again 

1 < e -^ 2 e ^i\\ 

for all (GC such that |ImC| > 7 and (u) ,9) 2 > 1 - /3 2 / 4 - This completes 
the proof of the Proposition. □ 

2.3. A quantitative Helgason-Holmgren theorem. Now we state and 
prove the main result in this section, the quantitative version of microlocal 
Helgason-Holmgren theorem: 

Theorem 2.5. Let M > 1 be constant. Given yo G R n , ujo G S n_1 , a > 
and j3 G (0, 1], consider (—a, a) xTcRx S n_1 introduced in (|2.1ip above, 
and define the dependence domain of the Radon transform data 

E = {2; G R n : (u,x — yo) = s, s G (—a, a), to G T}. 

Assume that for some p, 1 < p < 00, and X, < A < a function q 
satisfies the following conditions: 

(a) IeQ G XnL°°(R ra ) , where 1e stands for the characteristic function 
of the set E, furthermore 

\\<i\\l°°{e) + P-Eq\\x < M - 

(b) yo G suppg and suppg C {x G R n : (x — Uq,ujq) < 0}. 
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(c) (\,p,p)-Besov regularity on the dependence domain 
\\l E q- -y)Hi P(R n) 



r« \y\ n+Xp 



dy < M p . 



Then there exists a positive constant C = C(M, \G\, a, f3, A), such that 
(2.16) \\q\\ LP{G) < C 



log / (l + \s\r\\n yo q(s,-)\\ L1{r) ds 

J (—a, a) 



where 



(2.17) G = \ x e R n : \x - y \ < 

[ 8cosn(87r/p) 

A precise value of C in ()2. 16j) can be given as 

(2.18) C = C n M max (l, \G\* \ (1 + \y \) (a~ n + p~ n + a x 



Compared to Helgason support theorem, we relax the decay condition 
to the one given in ()2.14p and the L°° moduli of continuity are relaxed to 
integral moduli of continuity, which, under the condition < A < l/p, allow 
non continuous functions and are preserved (modulo constants) by multipli- 
cation by rough characteristic functions. These facts will be important in 
the applications. 

It will be convenient to use the classical 

= max(£, 0) t- = min(t, 0) 

to denote the positive and negative parts of a real number (or a function). 
The bounds on the Segal-Bargmann transform can be improved whenever 
the function / is supported on one side of a hyperplane. Indeed if 

supp/ C H£(s,ui ) 

then we have, for any yo S Hq(s,uo), 

(2.19) \%f(C)\ < (27Th)^e^ lm ^ 2 -^ Rc <- yo ^ 2 +\\f\\ L ^. 

Our first step will be to extrapolate estimate (|2.15p to capture points £ S C n 
with Im£ = 0. Note that for those values of the parameter the Segal- 
Bargmann transform is a gaussian transform. As in [5T] (see also \27\ Lemma 
9.6.5]) the following maximum principle for subharmonic functions will be 
the keystone of the proof. 

Lemma 2.6. Let a,b and A be positive constants. Consider 

R = {z G C : |Rez| < a, \Imz\ < b + e}, 
for some e > 0. Let F be a subharmonic function on R such that 

F{z) < (Rez_) 2 , 

for all z G R and 

F{z) < -A, 
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for z G R such that \Imz\ > b. Then, for 



|Im z\ < b, |Re z\ < — , 

we have 

F ( z ) < 

2a cosh (tt|) 

where 

A a \ 



(5 = min 



2a cosh (tt|) ' 3 



Proof. The claim follows by comparison of the subharmonic function F(z) — 
S 2 with the harmonic function 

e «--^|{j$ -»(;<* + «>)■ 

where 2 = x + iy, is in the rectangle Rg = [—5, a — 5] x [—b, b\. In fact, on 
the boundary of Rs we have 

F(x ± ib) - o" 2 < -A < G(ac ± ib) for x G [-8, a - 5], 

F(-S + iy) - 5 2 < = G(-<f + iy) 

and 

F(a - 6 + iy) - 5 2 < = G(a - 5 + iy). 
From the maximum principle F(z) < G(z) + <5 2 in Rs, which means that 

,2 x COsh (^) „,„ 



F(a + iy) < o" 2 - A -4f-f sin -(x + 5) 

cosh (^6) Va 

Since sint > 2t/ir for < i < vr/2, one has that 

9 cosh (^y) 2, 
F x + iy) < 5 2 - A— -|f|(-(x + 5), 
cosh (-6) a 

whenever < x + <5 < a/2. So if x is restricted to |x| < 5/2, then 

A 



F(x + iy) < 5 Z 



a cosh (S&j 

< 7 

~ 2a cosh 6) 

This completes the proof of the Lemma. □ 

Proposition 2.7. Consider q G XnL°°(R ri ) and Zei yo G R™ andwo G S"" 1 
be such that yo G suppg and suppg C {x G R™ : (x — yo,ujo) < 0}. Given 
a > and /? G (0, 1] consider the set 

r = { W £S n " 1 :(w, Wo ) 2 >l-^}. 
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// one has 

[ (l + \s\r\\n yo q(s r )\\ Ll{r) ds<e~^, 

J (—a, a) 

there exists a positive constant C, only depending on n, such that 
(2.20) e-^ Im <l 2 |T^(C)l 

< CM q (l + \ m \ + (T (1 + \s\) n \\n yo q{s, Oliver) ^ 



with 



K < — / \ n 2 ' M i '■= m ax(l, ||g|| LOO(Rn) + ||g|| x ), 

cosh f 



h 



a 2 



8|log/^(l + | S |)-||^ g( S ,-)IL 1(r) d S |' 
for all C G C n suc/i i/iat 

(2.21) |ReC-y |<^ r7^T> l Im CI < 



8cosh(8vr//3)' 1 (4-/3 2 ) l/2 ' 

Proof. Let £ G C n and denote z = {ojq,( — yo) G C. We write £ = (0 + 
(wo,?/o)Vo + u; with w G C n such that (Rero,wo) = (Imw,wo) = 0. Let us 
denote 

1= [ (l + \s\r\\n yo q(s,-)\\ LHr) da. 

J {—a, a) 

Choose 7 = ^ > as in Proposition \'2.'6\ then 
\%q{{z+ {uj ,y ))uj + w)\ 

< CM q {l + p+\y \) n e^ Imz ^ 2+ ^ lm ^ 2 (h~^l + h~^e-^ £ T 
< CM q {\ + p+ \y \) n e ^\ Imz \ 2 +^\ lmw \ 2 e~^4 

(_n J_o^. _n 1 3a 2 \ 

h 2 eah i6 X + h 2 e 2h is j 

for all z G C and 10 G C n such that 

(2.22) |Rez| 2 + |Re^-y + (^o,2/oH| 2 <« 2 /4, 

(2.23) |Imz|>2a//3, 

(2.24) |z + (w ,yo)| 2 + M 2 <p 2 , 
and 

(2.25) |Imz| 2 /(|Imz| 2 + |Imw;| 2 ) > 1 - /3 2 /4, 
where /? > is large enough. We next consider w G C n such that 

(2.26) \Rew-y + (coo,yo)u \ 2 < 3a 2 /16 
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and 

(2.27) \lmw\ 2 < 4a 2 /(4-/3 2 ). 
Then we have 

(2.28) \T h q((z+(uj ,y ))u + w)\<CM q (l+p+\y \) n h- n / 2 

x e ^(|i m ,p-|Re,| 2 ) +5 k|imH 2 x + e -ro^) , 

and conditions (|2.22p - (|2.25|l reduce to 

- 2 4a 2 . „ , 2 2 4a 2 



(2.29) |IM 2 <-, _<|T m ^< p 4 _^, 

with p large enough. 

Whenever X < e s , one can choose h = 8 | £ (0, 1] so that 

e2h 8 X = e 2h 8 . 

From (I2.19P we have 

(2.30) |7Jig((* + (wo,i«)»wo + u;)| < CM q h^e^ 



a _L|i mz |2 + i | Imt „|2_ i (Rez+)2 



for all u; G C n and zeC. Recall that M g = max(l, ||5||x,oorRn) + IMIx> 
Consider the sub-harmonic function $ defined as 



$(z) = \Rez\ 2 -\lmz\ 2 + 2h log \%q((z + (w , y ))w + u>) 

/ i 
e 2h i 



I 2 



+2/ilog 



C^/ 2 M g (l + p+|y |) n y ' 

where the variable w has been frozen. From (I2.30P and (|2.28[) one derives 
that 

$(z) < (Fiez_) 2 , 

for all z G C; and 



<&(z) < /ilogX 



a 2 



for all z G C satisfying ([2T29D . 

It is clear that $ satisfies the conditions of Lemma 12.61 with parameters 
a = b = j£ and A = hence we might conclude that 

hlog(J° a (l + \s\r \\n yo q(s, .)\\ Ll(r) ds)- 1 

<&(z) < — 

2(cosh(8vr//3)) 2 

for z G C such that 

a 2a 
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Choosing p = 4a/ (3 + 1 yo | , one can translate this estimate into the statement 
of the proposition, notice that (12361) . (|237j) and (pT29j) follow from (|23Tj) 
since 

|ReC - yol 2 = |Rez| 2 + \Rew-y + (ui ,y )u} \ 2 

and 

|Im£| 2 = |Imz| 2 + |Im«;| 2 . 
This completes the proof of the proposition. □ 



Proof of Theorem \2.5[ The key point is the fact that the Segal-Bargmann 
transform restricted to real values is a convolution with the Gaussian. We 
exploit this by means of the following lemma. Since actually we need a 
backward estimate for the heat equation, which is an ill posed problem, we 
requires the uniform Besov control of the potentials. 

Lemma 2.8. Consider q £ L p (R n ) and G an open set in R n . Assume that 
there exists A G (0, 1) such that 



Ilp(R" 



1/p 



\n+\p 



dy 



< +oo. 



x 'r- \y\ 

Then, there exists a positive constant C , only depending on n, such that 
\\q\\ L p(G) <c(hT% \\T h q\\ LP{G) + L q h* 

for all h £ (0,1]. 
Proof. Since 



q(x) 



(2nh) 

almost everywhere in G, 



1 ^ / n 1 
wT h q{x) + 



(27T) 



e 2 1 



'(q(x) - q(x - hzy))dy 



R" 



\LP(G) 



< 



1 



(2vr/i) 2 
1 

+ 



a IIT/i^Ir" II lp(g) 



(2vr)2 7r 



e~\\y\ 2 \q 



q(- - h*y)\ dy 



Lp(G) 



Minkowski's inequality ensures that there exists a positive constant C, only 
depending on n, such that 



i e-^l 2 
\q - q(- - hzy)\ „ dy 
(2ir) 2 



lp(G) 
< C p h^ 



\q - <?(• 



Ilp(R") 



R" 



\n+\p 



dy. 



This completes the proof of the Lemma. 



□ 
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Now we return to the proof of the Theorem. In the case 
f (l + \s\r\\TZ yo (q)(s r )\\ Ll{T) ds<e-^ 

holds, then Theorem 12. 51 is a consequence of Proposition 12.71 and Lemma [2,8l 
for the function IeQ- On the other hand, if 

I (l + | S |r||^ (g)( S , -)[lxi(r) ds>e-^, 

the conclusion of the statement of Theorem 12.51 is obvious. □ 

3. First application: illuminating Q from infinity (BU) 

We use the solutions of the Schrodinger equation in the maximal domain of 
the Laplace operator in Q. These solutions, satisfying the support condition, 
were constructed by Nachman and Street in [37] in the context of [33], but 
the construction in the case of illumination from infinity, [11] , is easier and 
follows the same steps. We will collect some estimates for these solutions 
which are (some of them implicitely) contained in [11] and [37] . 

Let H(£l; A) denote the elements of L 2 (Q) such that their weak Laplacean 
also belong to L 2 (£i). 

Lemma 3.1 (Bukhgeim and Uhlmann|llj). Assume that dQ. £ C 2 . Then 
the trace maps 

tr u = u\ d n 

and 

trm = d u u\ an , 

defined in C°° have an extension, again denoted as trj, j = 0, 1 which is 
continuous from H(Q;A) to the Sobolev space H~i~ l l 2 (dVL). 

If we assume in addition that trqu G H 3 / 2 , then u £ H 2 (Q) and 

\\ u \\m(n) + ll tr Hltfi/2(,9n) ^ c (\\ u \\H(n-,A) + \\ tv ou\\H 3 / 2 (dQ))i 
for some constant C > 0. 

The proof can be found in [11]. Let us remark that the definition of the 
extended trace maps is based on Green's formulae for smooth functions. 
Consider u £ H(£l; A), then on one hand, for u G H l / 2 (dVL), we have 

(3.1) troii(w) = / (uAv — Auv) dx 

Jn 

where v S H 2 (Q,) is the extension 

(3.2) v\ 9 n = 0,d u v\ d n = u. 
On the other hand, for uj € H 3 / 2 (d£l), we have 

(3.3) tritt(u;) = / (uAv — Auv) dx 

Jn 



tr u) - (tr u 



tri v 
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where v E i7 2 ($7) is the extension 

(3.4) v\ 9 q =U},d v v\ 9 Q = 0. 
The generalized Green's formula reads as follows. 
Lemma 3.2. For u E A) and v E H 2 (U), we have 

(3.5) / (A — q)uv dx = / u(A — q)v dx + ( tr^n 
Jn Jn ^ 

3.1. The Dirichlet-to-Neuman map. The next step is to define the Dirichlet- 

to-Neumann map associated to the Schrodinger equation — A+g. To achieve 

a definition in a extended domain that contains the traces of the solutions 

u E H(Q; A) of the equation (—A + q)u = 0, we will need the following 

lemma (see [37]). 

We will denote H(d£l) the range of the map 

tr : H(Q;A) -> H~ 1/2 (dn), 
and also consider the space of solutions of Schrodinger equation 
b q := {u E L 2 (Q) : (-A + q)u = 0} C H(Q; A). 

Then we have 

Lemma 3.3 (Nachman and Street [37J). If q E L°°(fJ) and is no£ a 
Dirichlet eigenvalue of — A + g in f/ien i/ie toace map 

tr : & 3 -»• ft(0fi) 

is one to one and onto. 

We consider the inverse maps 

P q = Uq 1 : H(dn) -> 6 g 
and define the norm of Tl{dVL) as 

IMk(an) : = II-Po^IIl 2 ^) 

Then we have: 

Lemma 3.4 (Nachman and Street [37]). The map tr : H(U; A) -)• 
is continuous and, under the hypothesis of the previous lemma, the map 

tr : b q -»• H(an) 

is a homeomorphism. 

One can define the Dirichlet-to-Neumann map 

(3.6) A„ : H(0fi) -> H~V 2 {dn) 
as the map 

A 9 (^)=tr 1 (P 9 (^)). 
To be more precise, by ([331) for p E H(dQ) and V E H 3 / 2 (dQ), we have 



(3.7) (\<pi/>) = f n ( P iV A v ~ qPq^v) dx, 
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where v is the extension in (|3.4p . It would be desirable to construct the 
Dirichlet-to-Neumann map as a selfdual operator, unfortunately this can 
not be achieved, instead we have: 

Lemma 3.5 (Nachman and Street [37]). Let qj , j = 1,2 be L°° potentials 
so that is not a Dirichlet eigenvalue of —A + qj in Q. Then A q2 — A qi 
extends to a continuous map %{dQ) — > %{dQ)* . 



These lemmas can be found in [37]. Let us remark that we have, if 
G H(dSl), 

(3-8) ((A g2 -A gi )M^) = I PM(qi-q2)P q M^- 



n 

This formula is the starting point of the recovery of values of the error in 
the interior of the domain. Notice that 

W^-qi ~ A <? 2 llw(an)^w(on)* > ~ ^■q2\\H 1 / 2 (dn)-^H~ 1 / 2 (dn)- 

3.2. The partial data map. Given N C S n_1 open (which could be very 
small) and £ G N, We assume F C dQ, to be a neighborhood of <9f2_(£) for 
any £ £ iV and .B a neighborhood of <9f2+(£) for any £ G iV. 

Given two potentials satisfying the hypothesis of Lemma 13. 5\ we will 
consider the difference of their partial data measurements as 

(3.9) || - A q2 \\* B ^F : = sup{ ^(A 91 - Ag 2 )ip B <pf) }, 

where sup is taken over the set 
{(<Pb,iPf) : \\<PB\\-H(dn) = \\^B\\n(dn) = 1,¥>b € £'{B) and G ^'(-F 1 )}- 

3.3. Fadeev's special solutions. We collect and remark results in |11] and 
|37j . concerning the existence and a priori bounds of iJ(f2, A) solutions of 
the Schrodinger equation adapted to the support requirements of the partial 
data. 

The modifications of [37] were done in the context of partial data consid- 
ered in [33] , but the construction can be adapted to the case of [H] , the only 
point is to write the operator conjugated with the exponential weight in the 
appropriate coordinates so that it is given as perturbations of the laplacean 
with terms of the complex operators d z and <9j. The final output of this 
construction is as follows. 

We consider q is as in Lemma 13.31 , £ and £ unit orthogonal vectors so 
that £ € N. We write x = (x±, X2,x") with respect to an orthonormal basis 
{ei, e n } so that e\ = £ and e2 = C an d x" G R n_2 . 

Theorem 3.6. For r > 1 sufficiently large and g G C°°(R n_2 ) ; and B 
and F as in the statement of Theorem I j.ffl there exists a unique solution 
w T G H(Q; A) of the equation ( — A + q)w T = in Q, such that trow T G 
%{d£l) n £'{B) and which can be written as 

w T (x) = e T{ t +ii: > x) (g(x") + R(t, x)) 



22 



where 



\\R(r, -)|U 2( n) < ck\\qg\\ LHn) + r x / 2 || Ag\\ L2{n) ). 



The same is true changing r by — r and B by F. 

The key ingredients in the proof are boundary Carleman estimates and 
some orthogonality properties of the reduced data solutions. 

Proposition 3.7. Let q £ L°°(J7), there exist tq > and C > such that 
for all u £ C°°(Q) , u\qq = and r > tq 

(3.10) Ct 2 [ |e- T ^n| 2 dx + r f {x,0\e~ Tx{x '° d u u\ 2 dA 

) (A - q)u\ 2 dx - t [ (x,i)\e- TX <d v u\ 2 dA. 



3.4. Stability from the Dirichlet-to-Neumann map to the Radon 
transform. We will use identity (|3.8p together with the special solutions of 
Theorem 13.61 to prove for q = {q\ — 92) 1^: 

Proposition 3.8. For any g £ C°°(R n ~ 2 ) there exists C > which only 
depends on f2 and the a priori bound of \\qj\\L°°, such that 



sup 



[£,£]- 



g{x") [ q(x" + t£ + s()dtdsdx" 

JR 2 



< C{T-^ 2 + e CT 



91 



A 



q2\\B->F)\\9\\H 2 (n), 



where [£, C] denotes the plane spanned by ^ and £. 
Corollary 3.9. // we consider the open set in S™" 1 

M = U ieN [$] ± , 
then we have in natural coordinates of the Radon transform 



(3.11) 



sup 



g{r)lZq{r, 77) dr 



R 



< C(t 



-1/2 



+ e CT ||A 9l 



A 



q2\\B^F)\\9\\H 2 {-R,)- 



Proof. With the notation of the Proposition, fix r\ £ [£,C] with |r/| = 1. 
We write x" £ [CjC]" 1 " as x " = x> " + rr / where x'" £ £, C] an d choose 
appropriate <? independent of x w , g(r) = g{rrf) in the Proposition to get the 
control of the Radon transform on any hyperplane, obtained by translation 
of hyperplanes in the pencil through the origin that contains £. We write 
x = rn + x'" + t£ + s£ with x'" £ [7/, £, £] , so that dx" = dx'"dr. Hence 



sup 



g(rn) 



R 



R"" 1 

CT I 



g(x w + rq + t^ + s() dt ds dx'" dr 



+ e CT ||A 



A 



<? 2 II j b^,f)II5 , IIh 2 (r) 



2.°, 



Translating the above in the natural variables (r, rj) of the Radon transform 
gives the Corollary. □ 

Proof of Proposition \3.8[ By Lemma 13.61 there exist v\ £ b qi such that 
supp tro^i C B of the form 

Vl = e - T{ ^' x) (l + R 1 {T,x)) 

and v 2 £ such that supptroi>2 C F and 

v 2 = e^-^ x \g(x") + R 2 (r,x)). 

We might write (I3.8p . by using these solutions as 

<^(A 92 - A gi )(tr fi) tr f2^) = J «i(gi - 52)^2 dx. 

We also have 

/(A 92 - A gi )(tr «i) tr v 2 \ < ||A gi - A q2 \\* B ^ F \\tr vi\\ n{dn )\\tr v 2 \\n(dn), 
which from Lemma 13.41 gives 



tr v 2 



< C\\A qi - A q2 \\* B _^ F \\vi\\ H (Q.^)\\v2\\H(n;A) 



(A q2 - A qi )(tr Q v{) 
Since 

IKIIff(n ; A) < C{\\q 2 \\ L °° + 1)IMIl2(Q) 

< C(||<fc|U~ + l)supe^ + ^> {\\g\\ LHu) + ||i?2(r, 011^(0)) 

< C(||ffilk» + l)e CT |l5ll^(n)- 
In a similar way, 

||^i||if(Q;A) < C(||gi||i«. + l)e Cr . 

Hence 

< Ce Cr ||A gi - A g2 ||B_ K F||s r ||.H" 2 (Q)) 



(A 92 - A gi )(tr < 



tr v 2 



with constants only depending on f2 and the a priori bound assumed on 
We have 

/ qviv 2 dx= / g(l + R 1 {r,x))(g(x") + R 2 (r,x))dx, 
Jq Jq 

hence 

/ Vtqv 2 dx > / qg(x")dx - Cr~ 1/2 ||5||^2 (n) 
Jj] Jq 

where C only depends on f2 and the a priori bound of ||<Zj||z,°°. 

Putting all together, we get 



/ qg(x")dx <Ce CT \\A qi -A q2 \\ B ^ F \\g\\ H 2 {n) +Cr^ 2 \ 
Jq 



9\\h 2 (Q)- 



This ends the proof of Proposition 13.8 



□ 
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Theorem 3.10 (Stability of Radon transform). We have the following es- 
timate on q = (gi — 92) In 



I l M n\Kq(r,r))\*dr) 



n + 3 
4 



2 

ra+3 



< + e -||A gi -A g2 ||^)STl 



Proof. The estimate can be obtained by interpolation of (|3.1ip and the fol- 
lowing estimate for the Radon transform 

/ / (l+r 2 )^ 1 |^g(r,7 ? )| 2 drda(7 ? ) < ||<z|| 2 2 . 



This estimate can be found in [38]. □ 
3.5. End of proof of Theorem 11.31 

Proof. Our aim is to use Theorem l2.5l together with estimate (|3.12p . We need 
the supporting condition (b) in Theorem 12.51 To achieve this condition let 
us take rj S M, we know by translation that there exists yo £ suppg such 
that the hyperplane H yo (0, rj), see (|2.ip . which contains yo satisfies condition 
(b). 

Now M is a neighbourhood of rj in the sphere and from the previous 
Theorem we can control the Radon transform for u) £ M and s£l. 

We can take a > large enough, so that G in (|2.17p contains supp q (the 
j3 depends on the size of M). Then 

(1 + \S - (u, C - yo)|) n ll^2/o(<?l - -)IIli(M) ds 



n + 3 



n+3 



< C(supp q,0) ^ ^JJlZq(r,r))\'dr j da( V ) j 
The choice r = \ log ||A 9l — A^Hj^pl in the estimate (|3.12p gives 

/ (1 + \s - (w,C - yo)\T pt yo {q\ - ^2)(s,-)|| L i(M) ds ^ C(supp<?,/3) 

J(-a,a) 

n-1 

x ((C| log ||A 91 - A ?2 ||^|-V2 + ( ||A 31 - A 92 ||^)V^ 
This together with (|2.16|) gives the desired estimate 



1 9i -92|liP(n) ^ C| log I log ||A c2l -Aq 2 \\* B ^ F \ 



for small norms of the difference of Dirichlet-to-Neumann maps. 

This ends the proof of Theorem 11.31 □ 
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4. Second Application: illuminating f2 from a point (KSU) 

To obtain stability in the case of [33] , we could have followed the same ap- 
proach, but in order to prove stability with the usual H 1 / 2 ^®) — > H~ 1 / 2 {dVl)- 
norm of the Dirichlet-to-Neumann map we will recall on the -ff 1 -solutions 
of the Schrodinger equation constructed by Chung (see |17j). in his work on 
the magnetic case with partial data. 

As above, we consider two potentials q\ and q 2 being in L°°(J7) and such 
that is not an eigenvalue of (-A + qj) : H^(9.) n H(U; A) — > L 2 (Q) for 
j G {1,2}. Let A q . denote the Dirichlet-to-Neumann map corresponding to 
the coefficient qj. 

Let (-|-) denote the duality between ii" 1 / 2 (c?r?) and H~ 1 ^ 2 (dfl) and recall 
that A Q] : i/ 1 / 2 (9$7) — > H-V 2 (dn) is defined as 

\ 3 fj 9/ := / ( Vtij , Vu) dx + / qjUjvdx 

for any fj,g G H l l 2 {d£l), where Uj G H 1 ^) is the weak solution of the 
Schrodinger equation (—A + qj)uj = in Q with Uj\gn = fj and v G H 1 ^) 
with v\qq = g. Since is not a Dirichlet eigenvalue, A qj is a well-defined 
bounded linear operator. 

An integration by parts shows that, for any Vj G H 1 ^) weak solution of 
the Schrodinger equation 

(4.1) (-A + g>j = 
in fi, one has 

(4.2) (^{A qi - A q2 )(vi\an) v 2 \dn^ = j (<?i - q2)viv 2 dx. 

Let yo be a point out of ch(f2), the convex hull of f2, and consider N a neigh- 
bourhood of uq also away from ch(il). Consider a hyperplane H separating 
./V and ch(O) and let H + denote the semi-space delimited by H and contain- 
ing 0. Consider R a positive constant such that, for all y G N, f2 C B(y, R) 
and set £ a subset of S n_1 with non-empty interior and 

S C {6 G S n_1 such that for all y G N, y + R6 £ H + and y - R9 <£ H + }. 

Let F and B denote open neighborhoods, for any y G N, of the faces 
dQ,_{y) and d£l + (y) respectively. Consider x £ C co (dQ) satisfying x '■ 
dtt — > [0, 1], suppx C F and xk = 1 with 9Q-(y) C F £ C F~ e C F (x and 
-F £ may depend on y). Obviously, we can write 

(( A ?i - A q 2 )( v i\dn) v 2 \an) = ((A ?1 - Ag 2 )(t>i|dn) (1 - xMan ) 

+ (( A gi - A g 2 )(' t; ilsn) X^lan) • 

Let w 2 G H 1 ^!) denote the weak solution of (—A + q 2 )w 2 = in Q with 
w 2\dQ. = u l|sn- This implies (^(fi — w 2 )\qq G if 1 / 2 (9f]) (^ stands for the 
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unit exterior normal vector on <90), since d£l is smooth enough (see [llj). 
Therefore, 



((A 9l - A q2 )(vi\ dU ) 



V2\dn 



< 



d v (vi - w 2 )(l - X)V2 dA 



+ ||(A 9l - Ag 2 )(^i|an)||^-i/2 (F) Wxv2\\ H y2 {F) 



(Hq (F) is the subspace of i? 1 / 2 (<9Sl) whose elements have their support in 
F, H~ 1 / 2 (F) is nothing but its dual -more about these spaces can be found 
in H5|). 

Assuming supp vi\qq C B, one has 



(A 91 - A q2 )(vi\ d n) v 2 \an 



+ IIA 



< / \d v (vi - ^2)11^2! dA 

'd(l\F e 



qi - Aq 2 l| ll^lll^ya^ llxllcWa+e^Q) I Ml #1/2 (,9^) ' 

where ||-|| = denotes the operator norm from Hq /2 (B) to H~ 1 / 2 {F) 

and e is any positive constant (C°' 1//2+,E (9Q) is the space of (1/2 + e)-H61der 
functions defined on the d£l -a close subset of R n ). 
For future references, 



(4.3) K^-A^Hen) 



v 2 \an) < / \d u (vi - w 2 )\\v 2 \dA 

JdQ\F E 

+ C||A gi - A q2 \\ \\vi\\ H l/2 {B) \\ V 2\\ H l/2 {dn) 



Since we assume that is not an eigenvalue of (—A + gA : Hq(£1) n 
H(Q;A) — ► L 2 {n)we have that d u ( Vl - w 2 )\ dn G L 2 {dfl). 

4.1. Controlling the Radon transform. Estimate (j4.3|) will be the start- 
ing point to control the 2-plane transform of q € L°°(R n ), for q = (qi—q 2 )ln- 
In order to obtain this information, we need to plug in (|4.3p special solutions 
of the Schrodinger equation (|4.ip . These special solutions are a generaliza- 
tion of the classical complex geometrical optic solutions (or Fadeev solutions) 
and were constructed by Kenig, Sjostrand and Uhlmann [33J and recently 
by Chung [17]. More precisely, from [33] (see Lemma 3.4 in [19]) we have, 

Proposition 4.1. There exists a solution v 2 of (—A + q 2 )v = of the form 

v 2 (t) = e ^ 2 +*)(a 2 +r 2 (r)), 

where t is a large parameter, ip 2 : (R n \ N) x N — > R defined by 

<P2(x,y) = - log |x- 

ip2 ■ x A x £ — >■ R defined by 

' x- y 



\x - y\ 
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with ft = n^gTv (-H+ n{i£ R n : |x — y\ < R}) and d S n-i the geodesic dis- 
tance on S n_1 associated to the Euclidean metric restricted to S n . Fur- 
thermore, a% : £1 x N x E — )• C defined as 

.,(,, ,,»)= pi, -,-..(, - , d;:;j;;::g| ) , 

with any g : S n_2 — )• C smooth, r^ir) £ iJ (O) and satisfies 
r\\r2{r)\\ LKQ) + r 1 / 2 ||r 2 (r)|| i2(an) + ||Vr 2 (r)|| 

< C (llfi'lL^S™- 2 ) + II^S n - 2 3llL2(S«- 2 )) 

where A S n-2 is i/ie Laplace- Beltrami operator on S n_2 /or i/ie canonical 
metric on S n_2 . .Bene C depends on H^IIjwq)) on 17, on i/ie distance of N 
to the hyperplane H . 

Let us denote 

Z(y) = {are0fi:<x-y,i/(aO)=O}, 

and assume that E is a compact subset of dQ^{y) \ Z{y) with E D <9Q \ 5. 
Let us recall Proposition 7.2 in [17j . 

Proposition 4.2. There exists a solution v\ of v\ of (—A + q{)v = wn- 
ishing on E and having the form 

Vl {r) = e T ^ 1+i ^\ ai + n(r)) - e r 'o, 

where r is the same as above, ip\ : (R n \ N) x iV — >■ R defined by 

<Pi(x,y) = -<P2(x,y), 

Tpl : $7 x N x E — >■ R defined by 

i>i(x,y,6) = -ip 2 (x,y,9), 

a\ : x N x E — >• C defined as 

_ n — 2 

aifoj/.tf) = (2|x-j/-0-(x-y)0|) 2 , 
additionally T\{t) G i7 1 (ri) and satisfies 

where C depends again on \\q\ H^oom) ■ Finally, I : QxNxT, — ?■ C is smooth 
and it satisfies Re / = y?i — A; u>it/i /c(x) ~ dist(cc, E) in G, a neighbourhood 
of E with non-empty interior in R n , and b : Q x N x E — > C is tance 
continuously differentiable in f2 and suppo C G. 

Let us remark recall the definition of the 2-plane transform. Given y £ 
R n , and 9 and r\ unitary orthogonal vectors, we denote 

Rq(y,9,n) := / q(y + t9 + rrj) dt dr. 
JRxR 



2N 



This is just the integral of q in the plane {y} + [9,rj\, where [6,rj\ denotes 
the plane spanned by 9 and rj. Notice that there is some redundancy on 
variables, since it is enough to define the above for y £ [6,rf\ L , see [40] , 
We assume that N contains the ball -6(0, a) and that ||A ?1 — A 92 || < 1. 

Theorem 4.3. The following estimate holds 



(4.4) sup \\Rq(y,e,r))\\ H -2 {Se) 



< 



ClV 1 / 4 Hall 1/2 +e cr IIA -A ll 1/4 



where Sg = S n 1 n #~ L and u;e consider the measure da, the volume form on 
Sg associated to the canonical metric on Sg. 

Proof. We next plug in (|4.2p the solutions in the above Propositions and 
bound by below the absolute value of the term in the right hand side of this 
identity: 



/ qv\V2dx 


> 


/ gaio^dx 


Jci 




Jci 



L 2 (nnG) 



+ ll a 2|| L 2 (n) lki|| L 2 (f7) + ||oi|| L 2 (n) IMI i2(n) + IM| L2(n) Ik2]|x2(n) ) • 



The last inequality, identity (|4.2p , the properties of the solutions v\ and V2 
and (|4.3p imply that 



(4.5) 



qaia 2 dx < \\q\\ l 2 (S n ~ 2 ) + \\^S n - 2 9\\L?(S n -' 2 )) 

+ / \d v (vi-W2)\\v2\dA + C\\k qi -k q2 \\e CT \\g\\ H i, sn -2 ) . 
Jdci\F F 



Mind that r" 1 / 2 comes from ||e _Tfc || i 2(^ nG ), dependences of C has not 
changed and c only depends on 0. 

We are next going to estimate the boundary integral term in (|4.5p . In 
order to do so, we are going to choose F e in such a way that dQ \ F e = {x 6 



dQ:(x-y)- u(x) > e} with y G N. Thus, 
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f \d v (v 1 -w 2 )\\v 2 \dA = [ e r ^|^(T;i-«; 2 )||a 2 +r 2 |d^ 
Jan\F e Jdn\F e 

<C[ [ e 2r ^\d u (v 1 -w 2 )\ 2 dA] 



dfl\F e 

X (llfl , Hfl' 1 (S n - 2 ) + II^S"- 2 fi f |ll / 2( S n-2) 



<C[- [ (u,x-y)e 2r ^\d u (v 1 -w 2 )\ 2 dA 
lan\F £ 



e 



X (llfllH^S"- 2 ) + II^S"- 2 fi f H J L2( S n-2)^ 



We now focus on the integral boundary term on the last inequality. Note 
that 

(4.6) (-A + q 2 )(vi -w 2 ) = -q\uvi 

with (vi — w 2 ) \qq = 0. We are going to use the Carleman estimate with 
boundary terms proved in [33] 

Proposition 4.4. Let q 2 £ L°°{Q), there exist tq > and C > such that 
for all u £ C°°(r2) , u\qq = and r > tq 

Ct 2 [ \e~ TLp2 u\ 2 dx + T [ {x-y,v(x))e~ Tip2 d u u\ 2 dA 
Jn Jdn + 

< [ \e~ Tlp2 (A-q)u\ 2 dx-T [ (x - y, v{x))\e~ TLp2 d v u\ 2 dA. 
Then we obtain, 

- f {v,(x-y))e 2 ^ 2 \d u (vi-w 2 )\ 2 dA) 1/2 
£ Jdn\F e J 

< C (V 1 / 2 UgllioofRn) + e cr \\du(vi - w 2 )\\ L 2 (dn _ {y))/ 

Here the constant depends additionally on e and the distance of N to ch(fi). 
We finally look at the L 2 -norm on the boundary. Let w 2 G i? 1 (SI) denote the 
solution of (14.1j) with j = 2 and the following boundary condition w 2 \ qq = 
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xMn -w 2 )\ d n £ H l / 2 {dtt). Then, one has 
\\d v {v x - w 2 )\\ L 2 (m _ (y)) 

/ d v (vi - W2)W2 dA I 
Jan J 

f V /2 
A(vi — W2)W2 dx + / (V(vi — W2), VW2) dx J 

n Jn J 

/ qviu>2dx+ / q2(v± - w 2 )w 2 + (V(vi - w 2 ), Vro 2 ) dx I 

\ 1/2 

qv±W2 dx I 



1/2 



<2 



Notice that the square of the last term appears in (|4.2|) if we change V2 by 
-u^. Thus, 



- w 2 )\\ 2 L 2 {d ^_ {y)) = ({A qi - A ?2 )(«i| 9Q ) w 



Taking into account that suppu^lao C -F one has 

\\d u (vi - W2)\\ L 2 {m _ {y) ) 

< [C \\A qi - A q J \\vi\\ H i {n) \\w 2 \\ H i {n) ) 

< (C ||A gi — A g2 || \\vi\\ H i in) \\d v (v 1 -w 2 )\\ H i/2 (dn ^ 1 



< 



C ||A 91 - A 92 || 1/2 e CT (|h - W2\\ L 2 [n) + ||Afa - «; 2 )|| i2(n) ) 



1/2 



In order to bound ||c?i/(^i — W2)||#i/2(gQ) in the last inequality, see [TT]. Since 
— u>2 is solution of the problem (|4.6p we can bound as follows 

\\d u (v 1 -W2)\\ L 2 im _ [y)) < C\\A qi - Ag 2 || 1/2 e CT ||g|| LO o (R n) , 
where the constant C depends on H^H^oo^™ with j £ {1, 2}. Summing up, 

J \d v (v 1 -W2)\\v 2 \dA < C (\\g\\ Hl(Sn - 2) + ||A S n-2 5 || L2(s , 



>dQ\F ( 

lL°°(R») TC \\ ly Ql ~ xy l2\\ ' ll ( /llL°°(R«) ' 



xl^lkll^^ + e-IIA^-A^H 1 / 2 
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Therefore, we finally get 



j^qaia 2 dx < — ||<?|| L °o( R n) (j|#|| L 2( S n-2) + ||A S n-25|| i 2( S T.-2) 

+ C ( t~ 1 ' 2 IMU,™ + e CT ||A 01 - K,\\ 1/2 " 



(4.7) 



L°°(R n )+ e H A gi A «j2ll II9Hl°°(R") ; 

S'll/f^s™- 2 ) + I|As"- 2 S'IIl 2 (S"- 2 ) 

+ C || Aqj — A g2 1| e CT Ibll/z^s™- 2 ) • 

In the following lines, we are going to relate the left hand side of the above 
estimate with the 2-plane transform of q. Note that 

/ q a\a2 dx = 2~^ n ~ 2 ' > / q(y + t9 + rr))g(r)) dt dr da(r]), 

Jn JrxRxSs 

where Sg = {rj G S n_1 : {rj, 0) = 0} and da is the volume form on Sg 
associated to euclidean metric restricted to Sg. Note that given y G N and 
£ E, one has 



<C ( -1/4 II ||l/2 +e cr|| A _ A || 1/4 

1 'Ia'IIh^s' 1 - 2 ) + I|As«- 2 5 , IIl 2 (s«- 2 ) 



for all y G iV and G E. Obviously, 

II^M,-)llH- 2 (5 e ) < C^-^lkll^^+e^ ||A 91 - A 92 || 1/4 
for all y G A and (9 G £. □ 

We would like to have an expression similar to (|4.4|) . with a norm of the 
Radon transform IZq in the left hand side, suitable to apply the stability 
result in section 2. This can be achieved in the three dimensional case, 
since the 2-plane transform Rq(y,9,rj), the integral of q in the plane P = 
{y} + i s a reparametrisation of the Radon transform. 

From now on, we will restrict our analysis to dimension n = 3. As we 
pointed out before, the variable y in (|4.4p is redundant, since the natural 
coordinates of the 2-plane transform are (y, P) where P = [0,rj] and y G P . 
The variable y in (|4.4|) can not be restricted to P-^, due to the 77- integral 
involved in the partial Sobolev norm. 

We have to make sense of the integral in (|4.4p in the appropriate coor- 
dinates of the Radon transform given by (s,£j) G R x S 2 . We can write 
poinwise Rq{y,Q,rj) = TZq(s,u), where s = (y,co) and oj = 9 x rj (vector 
product). To simplify the notation, we will assume that £ = £45 is bounded 
by the two planes parallel to H at distance 4<5 of the origin. 

We introduce geodesic polar coordinates on § 2 in the following way: given 
ujq G § 2 we can find #0 £ £<5 an d r/o so that they form an orthonormal frame. 
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We will denote 770 = ei, &o = ^2 and ojq = e%. Let 9 = #((/?) = cos(/?e2 + 
s\mpe\ with < 5. Let 77 run the geodesic Sg according to ip = dist(r), e%), 
we have rj = r)(<p, ip) = cos </? sin ^ei — sin(/9sin^e2+cos , i/'e3. We will perform 
the following calculations on the time zone 

r = T^S) = {Ti(<p,il>) : M < 5,0 < |^| < vr}. 

We will need to write Rq(y, 9, rj) in the standard coordinates TZq(s, uj(p, ip)), 
where uj(ip, ip) = 6((p) x 77(1^?, ip). For a fix <p the map 77 — > oj is just a rotation 
of angle tt/2 on the geodesic Sg, hence oj(ip,^) = rj(ip,^ — vr/2). As a map 
from r to r this rotation collapses the segments ip = ±tt/2 to the points 
±wo. We then restrict the time zone to 

(4.8) f = T^eM = { v (<p, 1>) : M < 5, if, € J}, 

where J = {5/8 < < tt/2 - 5/8} U {vr/2 + 5/8 < < vr - 5/8}. 
In these coordinates we have: 

Lemma 4.5. (Local estimate, n = 3) Let B(0,a) C N. Then 

r S / f \ 5/6 



(4.9) / / / \Rq(y,9,r,)\^d^) d^dy < 

JB(0,a) J-S \Jj J 



C ( -1/4 II ||l/2 +eC r|| A _ A HIM 173 

° I ' II "II L°°(R n ) H 91 92 II ) 

where C only depends on 5 and a . 

Proof. We will obtain (|4.9|) by interpolation of the following estimates. 
(4-10) / f \\Rq(y,9(<p), V (<p,-)) \\ H -2 (J) d<pdy< 



B(0,a) J-S 



c(t-^ 4 \\a\\ 1/2 +e cr IIA -A ll 1/4 
1 ' liy IIl°°(r») ~ e ll q2\\ 



and 



(4.H) / f [ (\d^R q ( yi 9(ip),n(ip,rb))\ 
JB{Q,a) J -8 J J 

+ \Rq(y,9(<p),ri{<p,il;))\)dil>d<pdy < C\\q\\ L 2. 

Estimate (|4.10p follows easily from (|4.4p . To prove estimate (|4.11f) we start 
with the derivative term, we change the order of integration and write y = 
s9 x 77 + y' with y' G [0, 77], 

/ \d^Rq(y, <%),?7(<£, V>)) \ dipdtpdy 

B{0,a) J -5 J J 

/S r ra r 
/ / / \d^Rq{soj + y\ 6(<p), 77(99, ip))\dy' dsdipdtp, 
-<5 J J J -a JBu 

where oj = u(ip, ip) and B^ = B(Q, a) n [0, 77]. 



Since y' does not change the 2-plane X-ray transform, we write the above 



as 



Ca 2 



8 J. J 



\d^Rq(su},9((p),T)((p,ip))\ dsdip dip, 



This expression can be written in terms of the Radon transform as 

/6 r rot 
\ \ \d^nq(s,u(ip,^))\dsd^dip 
-8 J J J -a 

\d^Kq{s, ri(ip, if> - tt/2))\ ds d^ dip. 



6 J. J 

To simplify the notation, following [38], we will use the homogeneous of 
degree —1 extension of TZ(s,rj) to 7] £ R 3 , to write 

dfi d df] 
d^TZq{s, r)((p, tp - vr/2)) = V v TZq{s, rj)- — = —K(xq)(s, rj) ■ —, 

then (|4.12[) can be bounded by 



(4.13) 



Co? 



8 J J 



d 

—K(xq)(s, rj(ip, ip - tt/2)) 



From the fact that on f one has 6/8 < 1 I 
that this can be bounded by 



dipdtp I 



dip dip ds. 
cosV'l < 1 — (5/8, we obtain 



< 



a 



2 ret 



K(xq)(s,rj) 



da(r]) ds. 



By using Cauchy-Schwarz inequality, this is majorized by 



< 



a 



5/2 



d 

—TZ(xq){s,7]) 



2 \ V2 

da{rf) ds 



5 1 / 2 

< C||ft(a;g)|| H i( RxS 2). 
Since in dimension n = 3, we have 

\\T^(xq)\\m(RxS 2 ) < C||(l + M)<7lU 2 , 

this fact, together with an easier estimate for the zero order term in (|4.1ip 
give, 

[ [ [ ^RqiyOx^ex^ldrjdedy^CdWxlqW^ + WqW^), 

JB(0,a) Jee(~8,8) JSe 

this ends the proof ()4.1ip and the Lemma. □ 
Now we return to the standard coordinates. 
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Corollary 4.6. (n=3) In the conditions of proposition assume that 
-B(0, a) C N. Then 



(4.14) / 1 [ \Kq(s,u)\dtr(u)ds 

J -a/2 JS 2 



< q( -1/4 II II 1/2 + e «-|IA A ||l/4 



where 1Z denotes the Radon transform. 

Proof. Let us take the open covering of S 2 given by {r a (5)} a eA, see (|4.8p . 
where A = {(ujo, 0q) S S 2 x : {uq, Qq) = 0}. By taking a finite subcovering 
and a partition of unity subordinated to it, we might reduce to prove the 
statement locally for the sets r a (<5) C S 2 . Then, for the local coordinates 
given in the lemma, ui = r](ip, ip — vr/2), we have 



pa/2 p 

(4.15) / / \Kq(s,u)\da(u)ds 
J -a/2 Jf„ 



q pa/2 p8 p 

< - / / / \Kq(s,ri(<p,il>))\d(pdil>ds. 

° J— a/2 J -8 J. J 

For any y G R 3 such that (y, rj((p, tp)) = s we have in terms of the two-plane 
transform 

Kq(s,T}((p,il>)) = Rq(y,6(<p),ri(<p,ip)). 

By taking y G 5(0, a/2), so that y = sr] + y' with y' G n^ , and denoting 
5 a = 5(0, a/2) n r/- 1 , we can write 

TZq(s, T]{y, ij>)) = \B a \~ l [ \Rq(srj + y', 9, V )\ dy' . 



B a 



a/2 JT 



/ / \Rq(srj + y',0,r])\dy > dipdipds 

J J J t3 a 



Inserting this in (I4.15h . we have 

ra/2 

\TZq(s, u))\ da(co) ds 

/a/2 r5 
/ 
-a/2 J-S ■ 

<C(a,5) [ f [ \R(y,e,r))\dipd<pdy. 

JB(p,a) J -8 J J 

The Corollary follows from the local estimate, together with Holder's in- 
equality. □ 

4.2. End of the proof of Theorem 11.41 We want to use the Corollary 
14.61 together with Theorem 12.51 Assume xo G P, the convex penumbra from 
N, and let y G N, so that (xq — y,v(xo)) = 0. To simplify notation we 
assume y = and take a > 0, such that 5(0, a) C N. 
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Assume xq E suppg (otherwise there is nothing to estimate) then the 
plane Ho(0, v(xq)), see (|2.ip . is a supporting plane of suppg, as required by 
(b) of Theorem 12.51 and, from the Corollary, we have estimate (|4.14p . We 
need, as required by Theorem 12.51 to use the Radon transform 1Z X0 with 
respect to the affine reference with center at the point xq E P. Consider 
P = il' where 5 = sup{|z — y\ : z E Q,y E N}. Then the set of planes 
{Hq(s,uj) : \s\ < a/2, ui E S 2 } contains the set {H Xo (s,u) : \s\ < a/8, to E 
T} where T = T(i/(x ),f3) is defined in (j2lTp . 

This allows us to write for I = (—a/8, a/8), 

(1 + \s - {uj,(- y )\) J \K Xo q(s,uj)\ da(u) ds 
< C I 1 I \Kq(s, u)\ da (u)ds 

J -a/2 JS 2 

< c(t~\ \\a\\ l/2 +e CT IIA -A ll 1/4 ^ 

Finally, the choice r = log ||A gi — A ?2 ||b_>.f|, together with Theorem 12.51 
gives the estimate in Theorem 11.41 for G a neighborhood of xq E P. The 
claim for G a neighborhood of P follows by standard arguments. 
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